Draft version April 5, 2011 

Preprint typeset using I^^'T^]X style emulateapj v. 11/10/09 



CONSTRAINING THE QUADRUPOLE MOMENT OF STELLAR-MASS BLACK-HOLE CANDIDATES WITH 

THE CONTINUUM FITTING METHOD 

COSIMO Bambi 

Institute for the Physics and Mathematics of the Universe, The University of Tokyo 
Kashiwa, Chiba 277-8583, Japan 

Enrico Barausse 

Department of Physics, University of Maryland 
College Park, Maryland 20742, United States 

(Dated: April 5, 2011) 
Draft version April 5, 2011 

^ ABSTRACT 

Black holes in General Relativity are known as Kerr black holes and are characterized solely by two 
parameters, the mass M and the spin J. All the higher multipole moments of the gravitational field 
are functions of these two parameters. For instance, the quadrupole moment is Q — —J'^/M, which 
^—^ implies that a measurement of Af , J, and Q for black hole candidates would allow one to test whether 

these objects are really black holes as described by General Relativity. While future gravitational- 
wave experiments will be able to test the Kerr nature of these objects with very high accuracy, in this 
^ paper we show that it is possible to put constraints on the quadrupole moment of stellar-mass black 

^ ^ hole candidates by using presently available X-ray data of the thermal spectrum of their accretion 

disk. 

Q Subject headings: accretion, accretion disks — black hole physics — general relativity — X-rays: 

binaries 
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1. INTRODUCTION 



The most general stationary and axisymmetric black-hole ( BH) solutio n of Einstein's equations in a four-dimensional 
asymptotically flat spacetime is given by the Kerr geometry ( |Kerr|1963 1 . Today there are at least two classes of astro- 
physical BH candida tes: stellar-mass objects in X-ray binary systems (mass M ~ 5 — 20 Mq) ( jCowley(T992; ,Remillard 
fc McClintock 2006 1 and super-massive objects at the center of most galaxies (Af ^ 10^ — 10^*^ M(?,) ( Kormendy fc 



Richst one"T995^. The existence of a third class of objects, intermediate-mass BHs with M - 10-^ - 10^ Mq ( |Millerl 
& Colbert 2004), is still cont roversial because their d etections are indirect and definitive dynamical measurements of 



their masses are still lacking ( [Miller fc Colbert | |2004| . 

All these objects are supposed to be Kerr Bhls because they cannot be explained otherwise without introducing new 
physics. In particular, stellar-mass BH cand idates in X-ray binary systems are too heavy to b e neutron or quark stars 
for any reasonable matter equation of state ( Rhoades fc RufRni|1974 Kalogera fc Baym|1996 |. Observations of stellar 
orbits around the super-massive BH candidate Sgr A* at t he center of the Galaxy show t hat this object is too massive 
compact, and old to be a cluster of non-luminous bodies (Maoz 



object made of sterile neutrinos, gravitinos or axinos supported by 



Rezzolla 


2007' 120081) 


Broderic 


c et al. 2009) 



1998]) or a fermion ball ("Schddcl ct al. '2002 | [i.e., a n 
aegeneracy p ressure ( Tsiklauri & Viollierlll998|}^ 



lottola 2004; Chircnti fc] 
] Broderick fc Narayan||200^ 



solution of Ge neral Relativ ity is still elusive 
put forward by 

spacetime. Such a generic spacetime can be used to describe the gravitational field around a central obj ect, whatever 
its natu re, an d its metric can be expressed in terms of the mass moments and current moments Se ( Geroch||1970 



A framework within which to test the Kerr BH hypothesis was first 
Ryan (1995 1997a b), who considered a general stationary, axisymmetric, asymptotically flat, vacuum 



Hansen||1974 1. Assuming reflection symmetry, the odd M- moments and even S'- moments are identically zero, so that 
the non-vanishing moments are the mass Mq = M, the mass quadrupole M2 = Q and the higher-order even terms 
Mi, Mq, . . as well as the angular momentum Si = J, the current octupole S3 and the higher-order odd terms 
S5, St, . . .. In the case of a Kerr BH, all the moments Mi and Se are locked to the mass and angular momentum by 
the following relation: 



Mi + iSi^ M [i 



J_ 
M 



(1) 



This is the celebrated "no- hair" theorem ( Carter|1971 Robinson|1975 Chrusciel fc Lopes Costa|2008 ): an (uncharged) 
stationary BH is uniquely characterized by its mass and spin angular momentum, 'i'heretore, a measurement of the 
mass, spin and higher moments (starting with the quadrupole moment Q) of BH candidates would permit testing 
Eq. (nl) and therefore the Kerr-nature of these objects. 

Ryan's idea was to use future gravitational- wave observations of extreme-mass ratio inspirals (EMRIs, i.e., systems 
consisting of a stellar-mass BH orbiting a super-massive BH in a galactic center) to perform this test. EMRIs will 
be a key source for the future-space based detector LISA: because the stellar-mass BH will orbit the super-massive 
BH ~ 10^ times during EISA's lifetime, slowing spiralling in due to the emission of energy and angular momentum 
via gravitational waves, even a small deviation from the Kerr geometry will build up an observable dephasing in 
the gravitational waveforms, thus allowing one to map the spacetime of super-massive BHs with very high accuracy. 
Ryan's spacetime mapping idea originated a whol e line of research aiming at using EIS A's o bservations of EMRIs to 
test the Kerr nature of super-massive BHs ([Collins fc Hughes 2004; Vigcland fc Hug hes|[2010[ Apostolatos et al. 20091 
Lukes-Gerakopoulos et al."2010^ 'Glampedakis fc Bab ak||2006( |Gair et al.||2008t jKesden et al.||2005i iB araussc et alj 
2007; Baraussc fc RezzoUa 2008; Barack fc Cutler 2007^ and even General Relativity itself ([Sopuerta fc Yunes 2009 



Barausse fc Sotiriou 200 8j). Anoth er independent (and complementary) test of the no-hair theorem with LISA uses BI ' 
quasi-normal modes (Berti et al. 2009). Because the frequencies of these modes encode the multipolar structure ([T| 



of the Kerr geometry, they c arfoe used to test consistency with the Kerr solution and to distinguish it from boson 



stars ( Berti fc Cardoso[[2006 ) or gravastars (Chirenti fc Rezzolla 



2007) 



Besides these tests based on gravitational waves, there are other proposals using electromagnetic radiation. Con- 
straints on__yifi_guadrupole moment of the compact companion of radio pulsars can be obtained with timing measure- 
ments (Wex fc Kopeikin 1999). Astrometric monitoring of stars orbiting at milliparsec distances frorn S gr A* may be 
used to test the no-hair theorem for the super-massive BH candidate at the center of the Galaxy ( Will|2 008; Merritt et 



aI] [20T0| . A very promising way to measure deviations from the Kerr metric is represented by fiature observa tions of 
the "shadow" of super-massive BH candidates t hrough very long baseline interferometry (VLBI) experiments (Bambi 
fc Yoshida||2010a| [Johannsen fc Psalt"is][2010a|b [). The study of quasi-periodic vari ability in BH s pectra may a' 
the geometry of the spacetime around BH candidates ([Johannsen fc Psaltis |2010c ). Remarkably, Psaltis & Joh 



so test 



Johannsen 



(2010d) also shows that the data for iron Ka emission lines in thin accretion disks can already constrain deviations 



from the Kerr geometry. Although these measur ements yield much less accurate constraints than what will be achieved 
with LISA, and can be subject to critiques (see Titarchuk et al.[ ([2009), who show that iron-Ka lines with the same 
features as those attributed to BH candidates are observed also around white dwarfs), these data are available now, 
which is not the case for all the o ther tests reviewed above. Ho wever, because of the controversial interpretation of the 
origin of these lines, and because Psaltis fc Johannsen (2010dl) flnds a degeneracy between the spacetime's quadrupole 
and spin (i.e., similar shapes for the iron-Ka lines can be obtained with a Kerr BH or with a non-Kerr object with spin 
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and quadrupole slightly shifted from the Kerr values), it is important to explore other techniques to test the no-hair 
theorem with present data. 

In this paper we propose using the continuum spectrum of BH candidates, which has b een shown to be potentially 



a promising too l to tell Kerr BHs from specific alternatives such as g ravastars (Harko ct al. 2009a), BHs in Chern- 
Simons gravity fHarko ct al. 2010b), BHs in Horava-Lifschi tz gravity p arko ct al. 2010a) or cert am classes of naked 
singularities ( jHarko et al.^2009b ; I'akahashi fc Ifarada 2010 Kovacs fc Harko_2010 ). While these attempts highlighted 
some important ditterences between the spectra of these objects and those of Kerr BHs, they relied on specific models 
for the BH candidate, and did not investigate whether presently available data allow one to break the degeneracy 
mentioned above between the parameters of these objects and those of a Kerr BH [i.e., whether present X-ray data 
can tell the spectrum of a non-Kerr object from that of a Kerr BH with arbitrary J and M). In this paper we address 
both issues, (i) by considering a very general model for the BH candidate [i.e., one which allows its quadrupole 
moment to slighlty deviate from the Kerr value, thus approximately describing a variety of almost-Kerr objects), (ii) 
by comparing our model to present X-ray data, although in a simplified way, and (in) by discussing the sources of 
systematic error that might affect the data and that must be properly understood before one can perform robust tests 
of the no-hair theorem. 

In the range 0.1 keV - 1 MeV, the generic spectrum of a stellar-mass BH candidate is characterized by three 
components, even if their relative intensities vary with the object and, for a given object, with time: i) a soft X-ray 
component (energies < 10 keV), ii) a hard power law X-ray component with an exponential cutoff (energies in the 
range 10 — 200 ke V, photon sp ectral index in the range 1 — 2.5), and Hi) a 7-ray component (energies > 300 keV). For 



a review, see e.g. Liang (1998). The soft X-ray component is commonly interpreted as the thermal spectrum of a thin 
disk, while the exact origin 01 the other two components is not so clear. 



Geometrically thin and optically thick accretion disks can be described by the Novikov-Thorne model (Novikov 



& Thorne 1973). They are expected when the accretion flow is radiatively efficient, which requires a luixiinosity 
L 0.3 LEdd-: where L^dd is the Eddington limit. The emission is blackbody-like. Assuming that the inner edge of 
the disk is at the innermost stable circular orbit (ISCOQ the disk luminosity of a Kerr BH is determined only by its 
mass, M, the mass accretion rate, M, a nd the spin parame ter, a ~ J/AP. This fact can thus be exploited to estimate 
the spin of stellar-mass BH candidates ( Zhang et al.||199T| . This is the continuu m fitting method and at present has 



been used to estimate the spin parameter of a few stellar-mass BH candidates (jMcClin tock et al.||20Tol p| Basically 



knowing the mass of the object, its distance from us, and the inclination angle of the disk, it is possible to fit the soft 
X-ray component of the source and deduce a and M. 

In this paper, we compute the thermal spectrum of a geometrically thin and optically thick accretion disk around 
a generic compact object. We use a subclass of Manko-Novikov spacetimes (Manko fc Novikov |199 2), which are 
stationary, axisymmetric, and asymptotically flat exact solutions of the vacuum Einstein equations. In addition to the 
mass and the spin of the massive object, here we have the anomalous quadrupole moment, q. The latter measures the 
deformation of the massive object with respect to a Kerr BH: when q > 0, the object is more oblate than a Kerr BH, 
when g < 0, it is more prolate, while, for g = 0, we recover the Kerr metric. The value of q determines the radius of 
the ISCO and changes the high frequency region of the spectrum of the disk. 

In general, this makes the spectrum of the disk almost degenerate in a and q. However, only in the Kerr case the 
radius of the ISCO goes to M as a approaches 1. For g ^ 0, even a small deviation from the Kerr metric makes the 
radius of the ISCO significantly larger than M. Since current X-ray observations suggest that there are objects with 
small ISCO radius, one can in principle obtain interesting constraints on the value of q. 

The purpose of this paper is therefore to present a preliminary investigation, showing that X-ray continuum spectra 
can potentially be used to constraint small quadrupole deviations away from the Kerr metric, once all the physical 
effects have been included in the model and all systematics have been understood. In particular, our computation 
of the disk's spectrum does not include the effect of light bending. This is a simplification of our model and there 
are no reasons for the light bending to be negligigle with respect to the other relativistic effects (Doppler boosting, 
gravitational redshift, and frame dragging). Another subtle issue is the computation of the spectral hardening factor 
(here not discussed), which is another weak point of our approach. Our study has to be taken as a preliminary 
investigation and significant work has still to be done before the continuum fitting method can be used to obtain 
reliable constraints on the Kerr geometry around stellar-mass BH candidates. 

The content of this paper is as follows. In Sec. [2] we review the basic properties of a geometrically thin and optically 
thick accretion disk and how to compute its spectrum. In Sec. [3] and H we present the results of our calculations, 
respectively for a Kerr BH and for a generic object with q ^ 0. In Sec.^ we show how current observations can be 
used to constrain g, while in Sec. |6]we discuss the possible systematic errors that could affect the continuum fitting 
technique and therefore out analysis. Lastly, in[7]we report our conclusions. The Manko-Novikov spacetime is reviewed 
in Appendix [A] and the properties of its ISCO are discussed in Appendix [B] 
Throughout the paper we use units in which Gn = c = 1, unless stated otherwise. 

2. THERMAL SPECTRUM OF A THIN DISK 



The standard model for a geometrically thin and optically thick accretion disk is due to Novikov and Thorne ( Novikov 
fc Thorne|1973| ). In a generic stationary, axisymmetric and asymptotically flat spacetime, one assumes that the disk is 



^ Such an assumption is supported e ither by observational facts ([ Steiner et al.||2010a [| and numerical simulations ([ Shafee et al.||2008| 
|Penna et al.|2010^ (but see |Noble et al.| | |2010| l). 

^ tor super-massive BHs, the disk temperature is lower (the effective temperature scales like Af""'^^) and this approach cannot be 
applied. 
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on the equatorial plane, that the disk's gas moves on nearly geodesic circular orbits, and that the radial heat transport 
is negligible compared to the energy radiated from the disk's surface. From the conservation laws for the rest-mass, 
angular momentum and en ergy, one can deduce three basic equations for the time-averaged radial structure of the 
disk ( Page fc Thorne|l974[) . These equations determine the radius-independent time-averaged mass accretion rate M, 



the time-averaged energy hux F{r) from the surface of the disk (as measured by an observer comoving with the disk's 
gas) and the time-average torque W^{r): 

M = -27rV^Su'' = const. (2) 
^ir) = -^^f{r). (3) 




WUr)^ ^ ^ - /(r). (4) 



Here S is the surface density, is the radial 4- velocity, G is the determinant of the near equatorial plane metric in 
cylindrical coordinates]^ i?, L^, and are respectively the conserved specific energy, the conserved z-component of 
the specific angular momentum, and the angular velocity d(j)/dt for equatorial circular geodesies, and /(r) is given by 

fir) = r [E - ni,) (drL,) dp , (5) 



(E-nL,) 



where rjn is the inner radius of the accretion disk and is assumed to be the ISCO radius. More details are given in 
Appendix |Bj 

Since the disk is in thermal equilibrium, the emission is blackbody-like and we can define an effective temperature 
T = T(r) from the relation F = crT'*, where a = 5.67 x 10^^ erg s— 1 cm^ is the Stefan-Boltzmann constant. 
Neglecting the effect of light bending, the equivalent isotropic luminosity is 



L(.).^cosz£ expK../(;T)]-l' 

where we have written explicitly the Planck constant h, the speed of light c, and the Boltzmann constant k. Here, 
i is the angle between the distant observer's line of sight and the direction orthogonal to the disk, and rout are 
respectively the inner and outer radius of the disk, while v is the radiation frequency in the local rest frame of the 
distant observer and Ve is the frequency in the local rest frame of the emitter. These two frequencies are related by 
the redshift factor 

"-^-"4, (7) 



where — (—1, 0, 0, 0) is the 4-velocity of the observer and = (u* , 0, 0, ilu* ) is the 4-velocity of the emitter. Using 
the normalization condition (7^,^?/^?/^ = —1, m* can be obtained to be 

ul - , ^ . (8) 

Because the t- and ^-component of a photon's canonical 4-momentum are conserved quantities in any stationary and 
axisymmetric spacetime, we can compute the quantity k^/kt at infinity. The result is k^/kt = rsin(/)sini and the 
redshift factor turns out to be 



\/~9tt - '^gt4,^ - gctxp^"^ 

1 -|- fir sin (h sin i 



(9) 



With g, we take into account the special and general relativistic effects of Doppler boost, gravitational redshift, and 
frame dragging. 

3. SPECTRUM IN KERR SPACETIMES 

To begin with, we calculate the thermal spectrum of a geometrically thin and optically thick accretion disk around a 
Kerr BH. Here we have four free parameters determining the luminosity (|6|: the mass of the BH M, the spin parameter 
a, the mass accretion rate M, and the inclination angle of the disk with respect to the distant observer i. However, 
usually M and i can be deduced from independent observations (see Section p] for an example). 

The role of the spin parameter is shown in Fig. ll where we assume M = 10 Mq, M = 10^^ g/s, and i = 45° . In 
the left panel, we present the radial profile of the effective temperature and, in the right panel, the observed spectrum 
vL(y). For a < 0, we mean that the disk is counterrotating. Since we assume that nn = '''isco' spin parameter 
determines the inner radius of the disk: as a increases, r^^ decreases and we find warmer matter at smaller radii. 



-G = v'o^grrS^, where a is the lapse function. In Kerr spacetime in Boyer-Lindquist coordinates, 
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Radial profile of the effective temperature (left panel) and spectrum uL{v) (right panel) of a thin accretion disk in Kerr 



spacetime for different value of the spin parameter a. Ifere we take the mass M 
inclination angle j = 45°. 



10 Mq , the mass accretion rate M = 10 g/s, and the 



1e+39 



1e+38 




1e+39 



1e+38 



1e+32 
1e+15 




1e+16 



1e+18 



1e+19 



1e+17 
v[Hz] 

Spectrum vL{v) of a thin accretion disk in Kerr spacetime with a 
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: 45°. Left panel: 



: 5, 10, 15 Mq and mass accretion rate M = 10 g/s. Right panel: mass Af = 10 Mq and mass accretion rate A/ = 0.5 X 10 



Fig. 2. 
mass M ■■ 

10^^ g/s, and 2 X 10^** g/s. 

At larger radii the effective temperature is essentially independent of the spin parameter. Therefore, a higher spin 
parameter moves the peak of vL(y) to higher frequency and to higher values. 

Changing the BH mass while keeping M constanl]^ has two effects. For larger masses, the effective temperature 
decreases (T c>c M^^/^) and therefore the peak of spectrum moves to lower frequency. At the same time, the size of 
the disk increases, increasing the total luminosity. In the left panel of Fig. [2] we show the cases M = 5, 10, 15 Mq for 

a = 0.9. The role of the mass accretion rate is shown in the right panel of Fig. [2] It is clear that a change in M only 
changes the effective temperature. 

The viewing angle i determines the effective disk surface seen by the distant observer and the correction due to the 
Doppler boosting (for i = 0° there is no Doppler boosting]]^ In Fig. [s] we show the cases i = 5°, 45° and 85° for a = 
(left panel) and a = 0.99 (right panel). 

Lastly, we show the effect of rin and Tout on the shape of the spectrum. So far we have adopted the standard 
assumption that the inner radius of the disk is at the ISCO and we have chosen the outer radius Tout = 10'^ M. 
However, if ri„ were larger than the radius of the ISCO, it would affect the high frequency part of the spectrum, 
mimicking a lower spin parameter, see Fig. |4] On the other hand, assuming a larger outer radius of the disk, the 
spectrum moves the cut-off at lower frequencies, with no changes at higher frequencies. 

4. SPECTRUM IN MANKO-NOVIKOV SPACETIMES 

Let us now consider the more general case in which the compact object is not a Kerr BH. The gravitational field 
around a generic compact body can be described by the Manko-Novikov metric, which is a stationary, axisymmetric 
and asymptotically flat exact solution of the vacuum Einstein equation and has an infinite number of free parameters. 
The structure of the spacetime presen ts strong similarities with the 5 — 2 Tomimatsu-Sato spacetime (Tomimatsu 



& Sato 



where t 



1972 Kodama fc Hikida||2003 ) . Here we restrict our attention to a subclass of the Manko-Novikov solution, 
le compact object is determined by its mass M, its spin parameter a — J /M'^, and the anomalous quadrupole 



While this assumption is useful to single out the effect of a change in M, if a BH accretes at the Eddington rate one has M oc M. 
^ We remind the reader that here we neglect the effect of light bending, whose contribution would also depend on the viewing angle. 
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moment q which regulates the deviations of the spacetime from the Kerr geometry: 

Q ^Kcrr 
M3 ' 



q = 



(10) 



where Q and Qkch = —a^M^ are respectively the quadrupole moment of the object and that of a Kerr BH. In this 
work, we consider only spin parameters \a\ < 1, which is the allowed range in the standard Manko-Novikov solutior^ 
As for the value of the anomalous quadrupole moment measuring the deviation from the Kerr metric: for q > the 
object is more oblate than a Kerr BH; for q < it is more prolate; for q = the Manko-Novikov metric reduces exactly 
to the Kerr metric. Thanks to this property, this subclass of the Manko- Novikov solu tion, for which wc give complete 
expressions in Appendix IaI is a perfect tool to set up a null experiment (Hughes 2006) to test the validity of the Kerr 
metric and of the no-hair theorem: any experiment pointing at a significantly non-zero value for q would imply that the 
compact object under consideration is not a BH as described by General Relativity. This us e of t he Manko-Novikov 
metric has been put forward in gr a vitational-wave astrophysics, n amely in Gair et al. (2008). See Collins & Hughes 
(2004); Vigeland & Hughes (2010); Glampedakis & Babak (20061 for other metrics wliich reduce exactly to the Kerr 
metric when the equivalent ot our anomalous quadrupole parameter q is set to 0, and which can therefore be used to 
perform null experiments testing the Kerr geometryFl 

In Figs. [5) |6) and [?} we show the radial profile of the thin accretion disk's effective temperature and the spectrum 
vL{v) for a few values of a and q. We still assume M = 10 M0, M = 10^* g/s, and i = 45°. For a given spin 
parameter, the value of q determines the radius of the ISCO - see also Appendix^ Since the temperature of the disk 
is higher at small radii, a non-zero q produces corrections in the high frequency region of the spectrum, while at low 

® In the Kerr case, \a\ < 1 is the condition for the existence of an event horizon. However, in the more general case of a compact object 
made of some kind of exotic rnat ter. the maximum value of \a\ may in principle be larger than 1 (Bambi & Frecsc 2009 Banibi et al. 20091 



|2010||Bambi & Yoshida|2010b|c| |. As in the case of the Tomimatsu-Sato family, also the Manko-Novikov solution can probably be extended 
to describe objects with \a\ > i (Manko & Moreno 1997). However, if s uch fast-rotating objects are very compact, they are most likely 



unstable, at least for sma ll values o t q, due to the ergorcgion instability ( Pani et al.|2010 i. 

The "bumpy BHs" of |Collins &: H ughes (2004); Vigeland & Hughes (2(^r seem related to the Manko- Novikov general solution becaus e 
they too allow an arbitrary multipolar structure for the non-Kerr spacetime. The "quasi-Kerr" metric of [G lampedakis fc Babak] | |2006l l, 
instead, is an (approximate) solution of the vacuum Einstein equations only for small values of the spacetime's spin, and allows only 
quadrupolar deformations. For these reasons, bumpy BHs or Manko-Novikov spacetimes are preferable options to test general deviations 
from the Kerr metric. 
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Fig. 5. — Radial profile of the effective temperature (left panels) and spectrum vL{v) (right panels) of a thin accretion disk in Manko- 
Novikov spacetime for spin parameter a = 0.99 (top panels) and a = 0.90 (bottom panels) and for different values of the anomalous 
quadrupolc moment q. Here we take the mass M = 10 Mq, the mass accretion rate M = 10^^ g/s, and the inclination angle i = 45°. 

frequencies there are no changes. The effect is quite small for slow-rotating objects or counterrotating disks, while it 
becomes relevant, and actually non-negligible, for fast-rotating bodies and corotating disks. There are two reasons for 
this: a small deviation from q 7^ produces a larger variation in the radius of the ISCO for higher spin parameters - 
see Fig. |9]in Appendix [B] ~ and, because the ISCO is closer to the compact object as a approaches 1, the spectrum 
of the disk is more sensitive to small deviations in the multipole moment expansion. This is the contrary of what 
happens in the Tomimatus-Sato spacetimes, where for \a\ — 1 all the solutions reduce to an extreme K err BH and 
thus deviations from the Kerr metric are more relevant for low spin parameters ( Bambi fc Yoshida||2010a ). 



5. OBSERVATIONAL CONSTRAINTS: THE CASE OF M33 X-7 

As an example of how accretion-disk thermal spectra can already put significant constraints on the deviation q of the 
quadrupole moment of BH candidates from that of a Kerr BH (c/. Eq. (10)), we consider the ca se of M33 X-7. Thi s 
object is an eclipsing X-ray binary consisting of a BH candidate accreting from a companion star ( Pietsch et al. 1 2006 ), 
and it s orbital parameters and its distance are measured with the highest accuracy among all known BH binaries ( jUrosz 
et al.|[2007J (see TablejlJ. In particular, the BH candidate's mass is measured to be M = 15.65 ± 1.45M0, while the 

disk's inclination is i =74.6° ± 1° and the distance is d = 840 ± 20 kpc ( |Orosz et aL]|2007|). 

The accurate knowledge of Af , i and d allows the continuum fitting method ( /hang et al.||1997 ) to extract reliable 
information on the spin of the BH candidate. By essentially fitting the Chandra and XMM-lNewton spectra of M33 X-7 
with a relativistic accretion disk model depending on the spin a, the Eddington ratio t = ^boi/iEdd (where Lboi and 
-^Ed d = 1-2572 X 10'^^(M/M (7)) erg/s are the bolometric and Eddington luminosities) and the hydrogen column density 
Nn 
{cf. 



Liu et al. 
i'able 1 of 



(|2008||201 [) me asur ed the spin to be a = O.84±O.O50 The Eddington ratio is instead £ = 0.0989±0.0073 
Liu et al. 1 2008 ) and Liu et al.'f2010')). We notice that the errors on a and £ include also the (propagated) 
effect of the uncertainties on M, d and i (Liu et al. 2008, 2010). 

Our simple disk model depends on three parameters, a, q and the Eddington ratio H — ^boiZ-^Edd- The latter can 
be rexpressed as i = M / M-Edd{o,, q), where we define the Eddington accretion rate as 



MEdd{a,q) 



L 



Edd 



ri{a,q)c'^ 



(11) 



r/ = 1 — -Ejgco (''isco ) being the efficiency of the conversion between rest-mass and electromagnetic energy ( Misner et 



* We warn the reader that this is the revised value reported in Liu et al. ( 2010 1, which corrects a bug in the analysis of Liu et al. 1 2008 I 
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Fig. 6. — As in Fig.js] for a = 0.50 (top panels) and a = 0.00 (bottom panels). 
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al.|1973 ). Ideally we would then have to fit the observed spectrum of M33 X-7 with this 3-parameter model. However, 
because of the difficulties and subtleties of analyzing the real Chandra and XMM- Newton spectra and because of the 
simplified nature of our disk model, we resorted to a simpler approach. While a thorough analysis of the real data 
will be needed to determine the precises constraints on q, our simplified treatment will show that such an analysis is 
definitely worth being done as it would permit ruling out entire regions of the (a, q) plane. 

In particular, instead of comparing our disk model with the raw data, we compare it to the spectrum of a thin 
disk with £* = 0.0989 and inclination i = 74.6° around a Kerr BH with spin a* = 0.84 and mass M = 15.65Af0 
(these are the values measured for M33 X-7). While meaningful and reliable constraints on the parameter a and 
q can only be obtained by fitting the original X-ray data, we use here this simplified approach because ours is a 
preliminary investigation and our results are only meant as a qualitative guide for future more rigorous studies. The 
spectrum is calculated with the standard Novikov-Thorne model reviewed in section |3] The observational errors on 
the "measured" spectrum are then mimicked by using the estimated final errors on the spin {Sa — 0.05) and Eddington 
ratio {6£ — 0.0073). Because the Eddington ratio regulates the bolometric luminosity {i.e., the normalization of the 
spectrum) one has L^''"{u,a*J*) < L^''"{iy,a* ± 5a, i* + Si) and L^''"{v,a*,£*) > L^''"{iy,a* ± 5a, i* - 51). It 
therefore makes sense to define the error as 



a(^) 



max(zyLK'="(t/, a* ± 5a, £* + Si)) - uim{iyL^''''^{iy, a* ± 5a, £* ~ 5i)) 



(12) 



To determine the values of a, q and £ giving the best fit, one would then have to minimize the reduced x^j which we 
define as 



Xrcd 



X2 

N 



1 



N ^ 

i=l 



(13) 



where the summation is performed over N sampling frequencies i^i and where L^'^ and L^°" are calculated as 
explained in Sections [4] and [3] respectively. To simplify the analysis, here we assume that the Eddigton ratio is fixed 
to the measured value**, and therefore seek to minimize 

'iy,L^^{v„a,q,£*) -iyiL^'="{iyi,a*,£*) 
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Fig. [5] for a = 0.50 (top panels) and a = 0.90 (bottom panels) in the case of counterrotating 
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Binary System 



M/Mq 



a 



Reference 



Shafee et al. 1 2006[ l 
Shafce et al. (2006) 
McClintock et al. (20"06| 
Davis et al. ( 2006 ) 
Liu et al. (2008 201o"l l 
Gou et al. (2009| 



4U 1543-47 
GRO J1655-40 
GRS 1915-1-105 
LMC X-3 
M33 X-7 
LMC X-1 
XTE J1550-564 



9.4 ± 1.0 
6.30 ±0.27 
14.0 ±4.4 
5-11 
15.65 ± 1.45 
10.91 ± 1.41 
9.10 ±0.61 



0.75 - 0.85 
0.65 - 0.75 
> 0.98 
< 0.26 
0.84 ± 0.05 

-1-0.05 
0.07 
4-0.20 
0.28 



0.92 
0.34 



Steiner et al. I 2010b 1 



TABLE 1 

Published spin measurements of stellar-mass BH candidates with the continuum fitting method. 



We stress that this simphfied approach makes sense because in principle the Eddington ratio can be determined 
independently, by integrating the luminosity over all frequencies. 



Since the spectrum of M33 X-7 is well fit with an accretion disk around a Kerr BH ( Liu et al.[2008 2010), one may 
wonder whether a fit with an additional parameter is statistically justified. As already stressed, our viewpoint is that 
X-ray continuum spectra can be used to constrain the value of the anomalous quadrupole moment, precisely because of 
the data's small error bars. This is, in other words, a null experiment, i.e. one in which we seek to measure a quantity 
which we expect to be zero. The measurement is therefore not one of the quantity itself, but rather a measurement of 
its "error bars" around its expected zero value. 

In Fig. [s] we show the contour plots of log]^Q(x^j,jj), as a function of a and q and as given by Eq. (14), in which 



we choose = 41 sampling frequencies Vi equally spaced, in logarithmic scale, from 10 ° to 1.5 x 1(F^ Hz. Also, 
to calculate and L^'^" we assume rout = lO^M, which allows us to significantly reduce the computational time 
needed to produce the spectra with respect to a larger outer radius. The regions of parameter space which are viable 
present logio(Xred) < 0. As expected, Xicd presents a minimum around 



0.84, q = (due to Eq. (14) 



but also the surrounding "valley" is in agreement with the data (see the right panel of 



with Xiod < 1 exists for q < —0.3 and a > 



ig-|8|). 

0.5, separated 



is exactly there) 

Moreover, a larger "valley" (featuring a central "basin" 
from the first one by a saddle. 
The physical interpretation of these two allowed "valleys" is quite straightforward. They stretch across the red line 
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Fig. 8. — l°Sio(Xrcd)' defined by Eq. l |l4[ l, for the comparison between the spectrum of M33 X-7 and that of a thin accretion disk 
in a Manko-Novikov spacetime with spin a and quadrupole parameter q (as defined in Eq. |10| ). Instead of fitting the original X-ray 
data, we use a simphfied model for the spectrum of M33 X-7 (see text for details), and therefore the constraints a and q only have a 
qualitative meaning. The viable regions are those with logio(Xred) ^' hence ruling out roughly half of the (a, q) plane. Our naive 
assumption | |12[ | for the errors probably overstimates the real observational uncertainties (see text for details). If the error | |12[ | were too 
large by a factor ~ 3.16 (10), X^cd ^^o'^l'i decrease by a factor 10 (100), effectively restricting the allowed (a,q) to the regions with 
l°gio(Xrcd) < ~1 (l°glo(Xrod) ^ ~2). The red line denotes the (a, q) for which the efficiency r;(a, q) equals that of M33 X-7 (see text for a 
physical interpretation). 

in Fig. [sj which corresponds the {a,q) for which T]{a,q) — rj{a ~ a*,q — 0), where 77(0 = CL*,q = 0) is the efficiency 
of the Kerr model that we use to mimick the data for M33 X-7[^ This fact is easy to understand. In our analysis we 
assume that the bolometric luminosity, Lboi = q)Mc^ is fixed and given by £* Lcdd = = a*, Q =^ 0)Af*c^, where 
M* denotes the accretion rate of M33 X-7. The accretion rate M is constrained to be close to AI* in order for the 
Manko-Novikov spectrum to reproduce that of M33 X-7 at low frequencies. This is because M basically regulates the 
slope of the spectrum at low frequencies: from Eq. (joj) one gets L{h') ^ T at small frequencies, but T cx M^/'' because 
of Eq. (|3| and the blackbody assumption. Therefore, if M ~ M* one obtains that it must be 77(0, q) ^ rj{a = a*,q = 0). 



We stress that we have determined these two allowed regions u nder the conservative a ssumption ( 12 1 for the error a 
In Eq. ( 12 ) we basically assumed that the errors determined by Liu et ah] ( |2008 2010[) for £ and a were uncorrelated, 
which could result in an estimate slightly larger than the real observational errors. I'his is hinted at also by Fig. 8j 
If one assumes q — (i e., if one adopts the Kerr-BH hypothesis) Fig. [s] shows that the allowed spins would be 
0.65 ^ a < 0.95, whereas Liu et al. (2008 20101 find a* = 0.84 ± 0.05. If our naive assumption overstimated the real 

2 



observational errors by a factor VlO « 
(0,(7) to the regions of Fig. ^ where logio,_,^^, 

0.79 5, a 5, 0.88, similar to the interval identified by Liu et al. ( 2008[ [2010[ ). Likewise, an observational error 10 times 



3.16 



Xred would dccreasc by a factor 10, effectively restricting the allowed 
(jj) < — 1. One can see that for q = 0, logio(x?ed) < ~1 would indeed give 



smaller than Eq. (12) would constrain the viable models to the regions with logio(x^(,(j) < —2. 

Even with our conservative assumption for the errors, however. Fig. [s] shows that more than half of the (a, q) p lane 
is ruled out. Nevertheless, if systematic errors for M33 X-7 were larger than assumed in Liu et al. (2008 2010), the 
constraints might be considerably weaker. For example, if the errors were \/lO ~ 3.16 larger than our assumption, 
only the region with logio(Xred) > 1 would be ruled out. We discuss possible sources of systematic error in the next 
section. We stress, however, that the presence of significant systematics would not only jeopardize our test of the 
no-hair theorem, but would represent a very serious problem also for the spin measurements with the continuum 
fitting technique, even if the Kerr-BH hypothesis is adopted. (This can be understood by looking at Fig. [s] for q = 0: 
as can be seen, the allowed interval for a grows rapidly if the error increases.) 

6. POSSIBLE SOURCES OF SYSTEMATIC ERRORS 

The continuum fitting method is a very promising technique for probing the space-time of stellar-mass BH candidates. 
Nevertheless, it is important to keep in mind that there are sources of of systematic errors that still need to be 
understood in order to obtain robust estimates of the spin parameter (if one assumes the Kerr-BH hypothesis) or 
constraints on the anomalous quadrupole moment with this method. 

The main source of uncertainty is the estimate of the hardening factor, sometimes called color factor, fcoi- Because 



^ The redline disappears for —0.29 < g < —0.05 because of the discontinous dependence of the ISCO on (a, q) (see AppendixrBl 



11 



in the inner part of the disk the temperature exceeds 10^ K, non-thermal processes are non-neghgible and the spectrum 
observed by a distant observer is not the blackbody-hke spectrum computed from the disk's effective temperature T. 
The hardening factor takes this effect into account, by replacing T with the color temperature Tcoi = fcoiT, and its 
typical values are in the range fcoi = 1-5 — 2.0. The computation of the hardening factor requires a reliable model of 
the disk atmosphere and its importance has bee n alrea dy stressed in 'Li ct al.' ( 20051). Significant progresses to address 
this issue have been done in Davis et al. (20051 and in Davis & Hubcny (2006 ). 

The continuum fitting technique also assumes that the spin of the compact object is perpendicular to the inner part 
of the accretion disk to within a few degrees. For stel lar-mass BH candid ates in X-ray binary systems, we expe ct this to 
be true, on the ba sis of binary population synthesis ([Fragos et al.|20l"0 1 . While the Bardeen-Petterson effect ( Bardeen 
&: Petterson|1975 1 may also be responsible for this efiect, for young objects the times cale ne cessary to align the central 



part of the disk turns out to b e too long. However, there are also observational data ( Maccarone,2002|) and theoretical 
arguments ( Fragile et al.|2001 ) s uggesting that tilted disks may be possible. Thi s assumption will be checked by future 
X-ray polarimetry observations ( Li et al.|2009 Schnittman &: Krolik|2009 2010), such as the GEMS mission scheduled 
for 2014. 

In our current analysis, we have also neglected the effect of light bending, because this is just a preliminary study 
to determine whether the continuum fitting method can conceivably be used to constrain deviations from the Kerr 
metric. While the effect of light bending must be taken into account in a complete analysis of the observational data, 
it has been quite commonly neg lected in similar preliminary studies appeared in the literature ( jHarko et al.||2009a|b| 
2010a|b[ [Kovacs fc Harko|[20To| . 



7. CONCLUSIONS 

If current astrophysical BH candidates are Kerr BHs, their spacetime should be completely specified by two param- 
eters, namely their mass M and spin J. This can be tested by measuring at least three multipole moments of the BH 
candidate. While there are several proposals to obtain such a measurement with future experiments, in this paper we 
have shown that current X-ray observations of stellar-mass BH candidates in binary systems can already be used to 
constrain possible deviations from the Kerr metric. 

We have computed the thermal spectrum of a geometrically thin and optically thick accretion disk around a compact 
object with mass M, spin parameter \a\ < 1 and arbitrary anomalous quadrupole moment q. For g = 0, we recover the 
Kerr metric. The exact value of q determines the inner radius of the disk, changing the high frequency region of the 
spectrum. The effect is small for low spin parameters or for counterrotating disk, but it becomes important for higher 
values of a. In general, the sole analysis of the disk spectrum cannot completely determine q, because the spectrum 
is degenerate in a and q and therefore one would need an independent measure of a. However, for very fast-rotating 
Kerr BHs the ISCO radius and therefore the disk's inner radius becomes very small. Any deviation from q = (i.e., 
from the Kerr solution) makes the inner radius grow quickly. Since current observations suggest that the inner radius 
of the accretion disk of some stellar-mass BH candidates is close to the gravitational radius Rg = GM/c^, one can 
constrain the anomalous quadrupole moment of these objects very efficiently. 

In this paper we have considered a specific ex ample, the stellar-mas s BH candidate M33 X-7, whose estimated spin is 
a — 0.84 ± 0.05 if one assumes it is a Kerr BH (Liu et al. 2008[ |2010[). Since stronger constraints on q can be obtained 
from objects with higher a, we could have c onsid ered GRB 1915 + 105, whose spin parameter has been estimated 
to be larger than 0.98 in McClintock et al. (2006) under the Kcrr-BH assumption. However, the measurements of 



the distance, mass and viewing angle of M33 X-7 are more reliable, thus making this object more suitable to obtain 
preliminary constraints on the a — q plane. 

To move our analysis beyond the simplified and preliminary stage we achieved in this paper, it is of paramount 
importance to properly understand the systematic errors that might affect the continuum fitting method, and which 
could in principle blur the difference between the spectra of Kerr BHs and those of other objects, and affect the 
measurements of the spin even if one adopts the Kerr BH hypothesis. Moreover, we will have to amend our disk model 
by including the following ingredients: 

1. The effect of light bending. A rigorous computation of the spectrum requires to trace the light rays from the 
surface of the accretion disk to the distant observer in the background metric. The effect of light bending is 
presumably no less important than the other relativistic effects and further alters the observed spectrum. 

2. The spectral hardening factor. In the inner part of the accretion disk, the temperature is high and non-blackbody 
effects cannot be neglected. We thus need an accurate model of the disk atmospher e for computing the spectral 
hardening factor ( Shimura fc Takahara||1995 Merloni et al.|[2000 Davis et al.||2005 ). 



3. Additional effects to be considered in an accurate study are the ones of limb darkening and of returning radiation. 



We would like to thank M. C. Miller and R. Takahashi for critically reading a preliminary version of this manuscript 
and providing useful feedback. E.B. would like to thank J. Brink and I. Mandel for pointing out a few typos in the 
original Manko-Novikov metric. C.B. wishes to acknowledge the Horace Hearne Institute for Theoretical Physics at 
Louisiana State University and the Michigan Center for Theoretical Physics at the University of Michigan, for support 
and hospitality when this work was finalized. The work of C.B. was supported by World Premier International 



12 



Research Center Initiative (WPI Initiative), MEXT, Japan, and by the JSPS Grant-in- Aid for Young Scientists (B) 
No. 22740147. E.B. acknowledges support from NSF Grants PIIY-0903631, and would like to acknowledge support 
and hospitality from the Institute for the Physics and Mathematics of the Universe at The University of Tokyo, where 
this work was started. 



APPENDIX 
MANKO-NOVIKOV SPACETIMES 



The Manko-Novi kov metric is a stationa ry, axisymmetric, and asymptotically flat exact solution of the vacuum 
Einstein equations ( Manko fc Novikov|l992 ). It is not a BH solutiorp"] but it can be used to describe the gravitational 
field outside a generic body like a compact star. The line element in quasi-cylindrical and prolate spheroidal coordinates 
is respectively 



=27 



-f (dt - ujd(i)y 
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1 and P„ are the Lcgendrc polynomials with argument xy/R: 
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(A2) 
(A3) 
(A4) 
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(A7) 
(A8) 
(A9) 

(AlO) 



(All) 



(A12) 



We notice that Eqs. ( |A6| ), (|A10|) and ( |All| cor rect a few typos in the original Manko-Novikov metric written in Manko] 
fc Novikov| ( [l992l ): see^rin^l^OOS^ a5d' |f^ ( [2007| ) . 



I'he solution has an inhnite nuniber ot tree pararneters: k, a, and a„ (n = 1, +cx)). For a ^ and q;„ — 0, it 
reduces to the Kerr metric. For a = a„ = 0, we find the Schwarzschild solution. For a = and a„ 7^ 0, we obtain 
the static Weyl metric. Without loss of generality, we can put ai = to bring the massive object to the origin of 
the coordina te system. I n this paper, we have restricted our attention to the subclass of Manko-Novikov spacetimes 
discussed in Gair et al. (20081, where a„ = for n 7^ 2. Therefore, we have three free parameters (fc, a, and 02), 



The Manko-Novikov spacetimes have naked singularities and closed time-like curves exterior to a horizon. Therefore, the no-hair 
theorem does not apply and the solution can have an infinite number of free parameters. Let us notice, however, that all these pathological 
features happen at very small radii and can be neglected in our study, because they are inside the inner radius of the disk. Here the basic 
idea is that naked singularities and closed time-like curves do not exist in reality because they are "covered" by some exotic object, whose 
exterior gravitational field is described by the Manko-Novikov metric. 
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which are related to the mass, M, the dimensionless spin parameter, a — J/M"^, and the dimensionless anomalous 
quadrupole moment, q ~ — (Q — QKm) / M'^ , of the object by the relations 

a = , k = Ai — — - , a2 = q-r^ ■ (A13) 

Let us notice that q measures the deviation from the quadrupole moment of a Kerr BH. In particular, since QKcrr — 
—a^M^, the solution is oblate for q > —a^ and prolate for q < —a^. However, when q 0, even all the higher order 
multipole moments of the spacetime have a different value from the Kerr ones. 

It is often useful to change coordinate system. The relation between the prolate spheroidal coordinates and the 
quasi-cylindrical coordinates is given by 



p = fcV(a;2 - 1) (1 - 2/2) , z^kxy, (AM) 

with inverse 



(A15) 



The relation between the standard Schwarzschild coordinates and the quasi-cylindrical coordinates is given by 

p = \/r^ - 2Mr + a^AP sin 61 , z = (r - A/) cos 6* . (A16) 
CIRCULAR ORBITS ON THE EQUATORIAL PLANE AND ISCO 
The line element of a generic stationary and axisymmetric spacetime can be written as 

ds^ = gudt"^ + 2gt^dtd(l) + QrrdP + gzzdz^ + g4>4>d4>'^ ■ (Bl) 

Since the metric is independent of the t and coordinates, we have the conserved specific energy at infinity, E, and 
the conserved z-component of the specific angular momentum at infinity, L^- This fact allows to write the t- and 
0-component of the 4-velocity of a test-particle as 

t Eg^^ + Lzgt4, ^ Egtcf, + L^gu 
u — — y , u = ^ . \^^) 

9U - 9ttgu 9U - 9ttg<t><t> 

From the conservation of the rest-mass, g^^u^W^ — we can write 

grrT^ + gzzZ^ = Kff z) , (B3) 

where the effective potential Vcs for fixed E and Lz is given by 

T. E'^g^^ + 2ELzgt4, + Llgtt , . 

VcS — T 1 • (o^) 



Writing the metric in quasi-cylindrical coordinates as in (All, one finds 

K=ff = ^ - 4 (Lz - LoEf - 1 . (B5) 

Circular orbits in the equatorial plane are located at the zeros and the turning points of the effective potential: 
r = i = 0, which implies VcS — 0, and r = z = 0, requiring respectively drVcs = and dzVes = 0. From these 
conditions, one can obtain the angular velocity, E, and Lz of the test-particle: 



d(j) -drgt4> ± \/ (drgt^f - {drgtt) idrg4,4,) 



^± = ^7 = 7^ , (B6) 

at Org^^ 

E = 9tt+9t^^ (37) 

\/-9tt - '2gt<t,^ ~ g^xp^'^ 



V-9tt - '^gt<i>^ - 300^^^ 

where the sign -f is for corotating orbits and the sign — for counterrotating ones. The orbits are stable under small 
perturbations if d^Vcs < and d^Vcs < 0. In Kerr spacetime, the second condition is always sa tisfied, so one can 
deduce the radius of the innermost stable circular orbit (ISCO) from d^VeS = . As first noticed in Gair et al. (2008), 
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Fig. 9. — Radial coordinate of the inner radius of the disk as a function of the anomalous quadrupole moment q for different values of 
the spin parameter a. Left panel: radial coordinate in the quasi-cylindrical coordinates. Right panel; radial coordinate in the standard 
Schwarzschild coordinates. Radial coordinate in unit M = \. 




-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 

q 

Fig. 10. — Radial coordinate in the standard Schwarzschild coordinates of the inner radius of the disk as a function of the anomalous 
quadrupole moment q for an extreme compact object a — 1. Radial coordinate in unit A/ = 1. 

in general that is not true in Manko-Novikov spacetimes. For q > 0, d^Vcs is always smaller than zero and one finds 
that the ISCO moves to larger radii as q increases. When g < 0, for any value of the spin parameter there are two 
critical values, say qi and q2 with qi > q2, and: 

1. For q > qi, one proceeds as in the case q > and finds the radius of the ISCO through the equation d^Vcff = 0. 
Tjgj^Q decreases as q decreases. 

2. For q2 < q < qi, there are two disconnected regions with stable orbits. The standard region r > ri, where ri is 
still given by d^Vcs = 0, and an internal region < r < r2, where r2 is once again given by d^Vcs = (that is, 
circular orbits with r2 < r < ri are radially unstable), while is given by d^V^s = (that is, circular orbits with 
r < are vertically unstable). As q decreases, ri decreases, while r2 and r^, increases. When q — q2, ri — r2. 



3. For q < q2, the ISCO is at r^, which is given by dlV^ 



off' 



0. 



For q > qi, the inner radius of the disk is the radius of the ISCO, i.e. rjn ~ t'isco- F'^^ 92 < q < qi, the energy and the 
angular momentum of the orbits in the region < r < r2 are higher than the ones at r > ri : the result is that the 
inner radius of the disk is at ri, since the accreting matter reaches the orbit at ri, which is a minimum of E and L^, 
and then plunges to the massive object. For q < q2, the inner radius of the disk is r^. 

In Fig. |9j we show the inner radius of the disk for a few values of the spin parameters a. In Fig. |10[ we show the 
case of a maximally rotating object with a 1: it is remarkable that only in the case q = the inner radius of the 
disk goes to M. For q 0, the inner radius of the disk is significantly larger. Since current observations suggest that, 
at least for some BH candidates, the inner radius of the disk is consistent with the one of a fast-rotating Kerr BH, 
deviations from the Kerr metric, if any, should not be large. 
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